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Abstract 

A detailed calculation of the real part of the finite temperature dynamic susceptibility of 
the free Bose gas is presented. After a short discussion on the different ways in which it can 
be calculated for temperatures above and below the Bose-Einstein transition temperature, 
its main properties and its evolution with q and T are analyzed. Finally, expressions for the 
lowest order energy weighted sum rules are also derived and studied. 
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The clear resolution of Bose-Einstein condensates in clouds of Alkali atoms at ultralow tem- 
peratures has revived the interest in the study of free and weakly interacting, homogeneous 
and inhomogeneous Bose systems during the last years Q. Recent experimental evidence ^] 
indicates that Bragg scattering can be succesfuUy applied to measure the dynamic structure 
function S{q,Lj) of trapped Bose-Einstein condensates, and hence to directly probe their mo- 
mentum distribution Q. 

During the last thirty years, much experimental and theoretical efforts have also been devoted 
towards the resolution of the condensate from neutron scattering in other, more correlated sys- 
tems like pure ^He and ^He-^He mixtures. However, none of the available formalisms have been 
able to accurately describe the influence of the temperature on the response, and particularly 
its influence on the condensate fraction, both at low and high momentum transfer 

In dilute systems, like the Bose-Einstein condensates of trapped alkali atoms, correlations 
are so weak that the response is expected to be qualitatively well described, at least at low 
momentum transfer, starting from the response of a free system and introducing the effect of 
the correlations through the Random Phase Approximation (RPA). In a previous work |^], a 
detailed description of the finite temperature dynamic structure function of the free Bose gas 
has been presented and discussed, focusing particularly on how q and T affect the total response 
and its coherent and incoherent parts. In this work we extend the analysis to the real part 
of the dynamic susceptibility xiQi^'iT), which is the building block from which the RPA is 
constructed. 

The dynamic susceptibility of an homogeneous free system is given by ^ 
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where p is the density, tp = /2m is the free kinetic energy spectrum, n{p) is the momentum 
distribution proportional to the occupation probability of each single-particle state of definite 
momentum, and v is the spin-isospin degeneracy which will be taken equal to unity throughout 
this letter. 

In the specific case of a free Bose gas case, n{p) is the sum of a condensate and a non- 
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condensate terms 
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where /? = 1/T is the inverse of the temperature and z is the fugacity, related to the chemical 
potential jj and the temperature through z = e^^ . The condensate fraction value no{T) varies 
with the temperature according to a T^/'^ law 
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where w 3.31p2/3/ 

m is the temperature at which the Bose-Einstein condensation sets in. 
Finally, the chemical potential /u(T) is fixed at each T by imposing the particle normalization 
condition 

J^Jdpnip) = l. (4) 

At finite temperature, the imaginary part of the dynamic susceptibility is related to the 
dynamic structure function S{q,u>;T) through the general relation 



Im [xiq, u; T)] = -vr [Siq, lo; T) - S{q, -co; T)] = -vr (l - e"^-) S{q, oj; T) , 
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where use has been made of the detailed balance condition, which relates the u; > and the 
< contributions to the finite temperature response in the form 



S{q,-u:-T)=e-P^S{q,u:-,T) . 
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The finite temperature S{q, uj; T) of the free Bose gas has been calculated in Q and yields, in 
terms of a new set of dimensionless variables that will be used throughout this work q = q/kB, 
oj = uj/eB and f = T/es where cb = k'^/2m and ks = p^^^ define the energy and momentum 
scales, 
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The first term in this expression takes into account the contribution of the atoms in the 
condensate and appears in the form of two delta peaks centered at the quasielastic recoil energies 
oj = ±q^. The second term describes the contribution coming from the rest of the atoms and 



shows itself as a function with two well defined peaks of finite width and height centered at the 
same energies. 

The real part of the dynamic susceptibility splits also in two pieces, corresponding to the 
condensate and non-condensate contributions, respectively. The condensate part Re[xc{Q, ^\ T)] 
is equal to nQ{T) times the real part of the T = dynamic susceptibility 
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thus diverging at lD = =bq^. On the other hand, the non-condensate part is a genuine contribution 
due to the finite temperature and is the main quantity analyzed in this work. 

The real part of the dynamic susceptibility can be calculated in different ways, depending 
on whether the temperature is above or below the Bose-Einstein transition temperature. In a 
standard treatment, one can just evaluate it by using a Kramers-Kronig relation 



Re 



-V I dCj' 



Im 



(9) 



UJ— UJ 

where V stands for Principal Value Integration. This expression relies on the implicit assumption 
that Im ti;;T) is an analytic function in the upper half of the complex plane, and that 
it decreases fast enough with increasing energy. Notice also that Eqs. (^) and (^) imply that 
Im[x{q,ijj',T)] and Re[xiq,i^',T)] are odd and even functions of w, respectively. 

Alternatively, one can evaluate Re x{q^Lb;T) as an infinite series, generalizing the result 
originally derived by Khanna and Glyde |9| for the free Fermi gas. Two different expressions 
are obtained in this way. In both cases, the dynamic susceptibility is written in the form 



x{q,Co;T) = x{q,0;T)+ I d^^x{q,tf) 
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where the partial derivative in the rhs can be evaluated from the definition in Eq. (|l]) and is 
expressed as an integral that may be carried out on the complex plane. As a result one gets 
an infinite series corresponding to a sum of residues, that has to be inserted in Eq. (|lO| ) and 
integrated to get the total dynamic susceptibility. 

In a first approach, one integrates this series term by term to end up with 
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where the index n runs from to cxd. In Eq. a„ and 6„ are the real and imaginary parts of 
the poles of the momentum distribution laying in the first quadrant 
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while yi and §2 are the usual West scaling variables, related to the momentum and the energy 
transfer through 

& m = l{^ + ^) ■ (13) 
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Alternatively, one can first evaluate the sum of all the residues contributing to the energy 
derivative of the real part of the susceptibility 
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and afterwards perform the integration. Although this method seems to be more elaborated 
than the simple sum in Eq. ([TT|), it is usually preferred since the former expression is known to 
converge rather slowly. In both cases, however. Re x{q, 0; T) is a required quantity that can 
be obtained, proceeding as before, in the form of a series 
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and the sum extending from to 00 with a„ and bn given in Eq. (12). 

Finally, Re[x{q,ob;T)] can also be evaluated by direct numerical integration of Eq. (|l]). After 
simple manipulations, one finds 
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which can be safely carried out since their singularities are of the integrable type. 

For temperatures above Tc, all three methods are equivalent and any of them can be used to 
obtain the non-condensate contribution to the real part of the dynamic susceptibility. Below Tc, 
however, only the direct integration method and the Kramers-Kronig relation can be applied, 
even though in the latter case the analyticity condition required for Im[x(q,uj;T)] is clearly 
violated since at those temperatures the chemical potential vanishes and the fugacity z goes to 
1, therefore the function presents two singularities on the real axis located at the quasielastic 
energies u) = iq-^. However, these are logarithmic singularities that can be integrated on the 
complex plane and yield no additional contribution to the integrals. On the other hand, none of 
the series methods mentioned above can be used because below Tc the momentum distribution 
grows as p^"^ for small p, thus making i?e[x(0, 0; T)] in Eq. ( [l6| ) diverge. 

Since jl(T < Tc) = 0, the temperature dependence of Re[x{q,(^;T < Tc)] can be easily 
extracted once it is written in terms of two new variables Q = q/Vt and v = uj/f. Moreover, 
one can define a new function 

Re[x{Q, i^)] = -^Re \x{q, f)] (18) 



'T 

that does not depend on the temperature when Q and v are taken as the new independent 
variables. Actually, this scaling property is also satisfied by the imaginary part of xiQ^^^T) 
which inherits it from the T <Tc behaviour of the dynamic structure function of the free Bose 
gas p. 

Re[x{Q,i^)] is shown in Fig. for several values of Q. At low u, Re[x{Q,i^)] is nega- 
tive and discontinuous at u = Q"^, jumping to positive values and afterwards decaying to 0. 
The discontinuity at = is a direct consequence of the divergent behaviour of the mo- 
mentum distribution at low momenta, and is therefore directly related to the presence of a 
Bose condensate. On the other hand, the decay of Re[x{Q,i^)] at large energies is of the form 
{I — nQ[T))2Q'^ / u'^T^/'^ = 2Q'^ /u'^Tc^'^ . This last property is mirrored from the fact that the real 
part of the dynamic susceptibility is related to the dynamic structure function through Eqs. (P) 
and (^, and this together with the /-sum rule satisfied by S{q,uj;T) is enough to determine 
the analytical behaviour of the real part of the susceptibility at large energies Q. Finally, the 
factor 1 — no(r) reflects the fact that we are only looking at the non-condensate contribution 
to the susceptibility, while the condensate term adds the extra strength required to have a total 
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decay of the form jUp' as expected. 

With rising momentum transfer the discontinuity shifts to higher Q's, the strength at low 
energies and around is depressed, but the energy range where i?e[x((5, i^)] presents a signifi- 
cant contribution is increased. Due to the scaling in Q and the effect of fixing T and rising 
q is equivalent to the effect of fixing q and decreasing T. Therefore, the evolution with T of 
i2e[x(Q', 'i; r)] at fixed q can be read in the figure by moving from higher to lower Q's, while the 
evolution with q at fixed T is represented by the same sequence of curves but in reverse order. 
At T = the non-condensate contribution to the dynamic susceptibility vanishes and only the 
condensate term is left. 

When T crosses the Bose-Einstein transition temperature, the discontinuity at lD = is 
smeared out and Re\x{f[^^'iTy\ becomes a continuous function of u), as shown in Fig. (|2|) for a 
momentum transfer oi q = \. This is due to the fact that, above Tc, the chemical potential is 
no longer and thus the p ^ divergence in the momentum distribution is removed. Actually 
the process is quite fast, and already at T = 7 and a momentum transfer of g = 1, no apparent 
trace of a discontinuity in the real part of the dynamic susceptibility is left. 

The evolution with T of the non-condensate contribution to the real part of the dynamic 
susceptibility for T > and for three values of the momentum transfer is sketched in Fig. (^) . 
The actual function depicted is Tiie[x(Q', 5"^; T")] with Y = {v/Q — Q)/2, since in the high tem- 
perature limit the susceptibility of the free Bose gas approaches the T ^ oo classical prediction 
which for energies around the quasielastic peak u = q^ can be casted in the form 

TRe[xci{q. Y;T ^ oo)] = -—V / dz ^ , (19) 

\/TT J-od z-Y 

thus being a function of Y alone. Since Re[x{q,oj;T)] is and even function of lD, only the 
contribution at positive energies is shown, so the initial point in each curve is 1" = —Q/2. 
As it can be seen from the figure, the departure from the classical limit at low temperatures 
is sizeable, while the Bose prediction gets closer to the classical one when the temperature is 
risen, as expected. However, the way in which the classical limit is reached depends on the 
momentum transfer. At low q, the low temperature susceptibility rapidly varies from negative 
values below the classical prediction to positive ones above it. With increasing momentum 
transfer, the strength of the Bose function at fixed temperature is progressively reduced, but 
the range of maximum variation in which the function goes from negative to positive values 
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grows. In all cases, however, the large Y tails of the Bose function conform to the classical 
prediction at any temperature, a fact that should not surprise since as commented above the 
dynamic susceptibility is known to decrease as Ic^ juP' independently of the temperature and 
the statistics. 

The behaviour of the real part of the dynamic susceptibility at fixed T larger than Tf. and 
for a sequence of increasing values of the momentum transfer q is depicted in Fig. In this 
case, the g — > oo limit of Re\)({(i^Y \T')\ at finite Y can be obtained from the leading term in a 
1/g expansion, a property that is inherited from the high q behaviour of the dynamic structure 



function [10|. As in the latter case and for the free Bose gas Q, the function 

lRe\x{q.Y-f)\ (20) 

approaches a fixed curve when g ^ oo that depends on the temperature only through the 
chemical potential. This limiting function is represented with a solid line in the figure, while 
the other curves show the way in which this limit is approached with increasing q. Notice 
that, contrary to what happened in the previous case, increasing q at fixed T makes the initial 
point Y = —Q/2 in the limiting curve go to — oo, thus emphasizing that the region of maximal 
variation of Re[xiq,Y;T)] lies always around y ~ corresponding to energies close to the 
quasielastic peak uj = (f' . Notice also that the limiting curve is always reached from the right, 
as {q/T)Re[x{q,Y;T)] shifts to lower values of Y when q is risen. 

The dependence of the real part of the susceptibility on T and q can also be analyzed by 
looking at the sum rules it satisfies, which are defined as the energy-weighted integrals 

rhn{q;f)= dCoCj^Re[x{q,Co-f)] . (21) 



Sum rules could also be defined extending the integration range to a) G (— oo, oo), but then all 
the odd order ones would trivially vanish because Re[x{q,^]T)] is an even function of a). With 
the previous definition, and taking into account that the real part of the dynamic susceptibility 
decreases as 2c^ jiiJ^ at large energies, one readily notices that no sum rule with n > exists. 
However, an expansion around u;~^ = displays a next-to-leading term of order cD~^, so one 
can still find a second energy weighted sum rule by subtracting the leading contribution at high 
energies. Proceeding in this way one finds the following relations 

mo(g;T) = / dLDi?e[x(g,cD;r)] =0 , (22) 
JO 



m2{q; f) = du (nelxiq, u; T)]-^] = . (23) 



These results actually apply to the total real part of the dynamic susceptibility, including 
the condensate contribution. However, they are general and therefore hold also for the T = 
susceptibility of the free Bose gas also. The condensate contribution at finite temperature is 
equal to no{T)Re[x{q,uj;T = 0)], and thus its 'mo{q;T) and rh2{q;T) also vanish. In this way, 
therefore, the above sum rules apply separately to the condensate and non-condensate parts 
of the dynamic susceptibility, respectively. Additionally, the fact that both sum rules vanish 
indicates that Re[xiq, T)] must have different regions of positive and negative values, a feature 
that is in fact general to all systems since the previous sum rules apply to both correlated and 
uncorrected systems at T = and finite temperature. 

Although the previous results hold for any system conserving the total number of particles, 
it is difficult to obtain generic expressions for higher order sum rules since then additional terms 

in the tZ-'^^ expansion of i?e[x(g, cD; T)] must be subtracted, and the coefficients of this expansion 
are related to sum rules of the dynamic structure function S{q,oj;T), which arc specific to each 
system. Furthermore, the simple subtraction of 2g^/a;^ does not allow for the derivation of a 
mi{q;T) sum rule, a problem that remains even when higher orders in the u;~^ series expansion 
oi Re[x{q,oj;T)] are subtracted. 

An alternative set of sum rules, valid for correlated Bose systems at T = and T > and 
for the noninteracting Bose gas at T > 0, can be derived by subtracting to xiQi'^i T) the T = 
dynamic susceptibility x^{q,uj;Qi) of the free Bose gas. As long as the real part is concerned, 
these become 

M„(g,T) = d^(^'' {Rem^u-f)] - Re\x\q,uM) ■ (24) 

Using the explicit representation of the dynamic susceptibility in terms of the dynamic 
structure function, the first M„ (g, T) moments can be easily derived and yield 



Mo(^;r) = 0, (25) 
Mi(g;f) = - j^^du S{q,u-f)\J'^ , (26) 



M2{q;f) = 0, (27) 
Interestingly enough, the expression of Mi{q; T) indicates that even though it is not know a 
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priori whether the dynamic structure function S{q, a); T) of a given system does have sum rules 



to all orders |11], still some combinations of them may exist. Notice also that this sum rule is 
trivially satisfied by the condensate term, since the condensate contribution to S{q, a); T) is a 
sum of two delta peaks centered at cD = itq^ where the logarithm cancels. Therefore and for the 
free Bose gas, the non-trivial part of Mi{q;T) must be satisfied by the non-condensate term 
alone. Moreover and from the same argument given above, it is also apparent that the total set 
of sum rules must be separately satisfied by the condensate and non-condensate contributions 
to Re[x{q,uj;T)]. Finally, the zeroth and second order sum rules Mo(g;T) and M2{q;T) are a 
direct consequence of the results shown in Eqs. (^) and (pSj). 

The fact that below Tc the temperature dependence of both the dynamic susceptibility and 
the response function can be extracted as indicated in Eq. (^) implies that an energy weighted 
sum rule can also be defined in the form 

M,iQ;T) = ^M,{q;T)^ J_^d,^ S{Q,u;T)\n i^j , (28) 

where S{Q;i';f) = Snc{Q,'^',f)/\/T is the non-condensate contribution to the response of the 
free Bose gas written in the variables Q and i/. Below the Bose-Einstein transition temperature, 
this factorization removes any dependence on T and thus Mi becomes a function of Q alone. 
Above Tc, however, a dependence on T remains through the chemical potential /i(T'). Mi(Q;T) 
is shown in Fig. (^) for several temperatures below and above Tc- In all cases, the low Q 
behaviour is divergent, a feature produced by the term inside the logarithm. On the other 
hand, the high Q behaviour is dominated by a slowly decreasing tail that approach from below. 
In the intermediate region, Mi(Q;r) presents a minimum at some Q that increases with the 
temperature, even though it is hardly noticeable at high T's. 

In summary, it has been shown that the non-condensate contribution to the real part of 
the dynamic susceptibility of the free Bose gas at finite temperature can be evaluated in three 
different ways: by direct integration, using a Kramers-Kronig relation and performing a Khanna 
and Clyde -like series expansion. While the first two methods can be used above and below the 
Bose-Einstein transition temperature, the latter fails below Tc since at those temperatures the 
starting point in the series diverges. 

Rs[x{Qt'^tT ^ ^c)] is shown to scale in two new variables Q = (j/Vt and u = uj/T from 
where its temperature dependence can be extracted, a condition that is lost above Tc due to the 
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finite value taken by the chemical potential fi{T). Below Tc, the real part of the non-condensate 
contribution to the dynamic susceptibility presents a discontinuity at a; = ibg^ produced by the 
singular behaviour of the momentum distribution at p ^ 0. When T exceeds Tc this singularity 
in n{p — 0) is removed and as a result the previous discontinuity is smeared out. At high 
T > Tc, Re[x{q,oj;T)] approaches the classical prediction computed from a Maxwell-Boltzmann 
momentum distribution which for energies close to the quasielastic peak {j = can be brought 
to a form that scales in a single variable Y = [tj/q — q)/2\/T. On the other hand, the leading 
contribution to the high q limit of Re[x{q,oj',T)] can also be expressed in terms of Y alone, a 
property that is inherited from the large momentum transfer behaviour of S{q,Lb;T). 

Low order energy weighted sum rules of the non-condensate contribution to the real part 
of the dynamic susceptibility are also derived and discussed. Since the long energy tails of 
Re[x{q,oj;T)] decrease as l/Co^, no sum rule of order higher than exists. However, a new set 
of sum rules, obtained by subtracting the T = dynamic susceptibility of the free Bose gas, is 
proposed and analyzed. These sum rules are general and apply to any system conserving the total 
number of particles, and their analysis show that in general the real part of the susceptibility 
must have different regions where it changes sign. Despite the simplicity of the model analyzed, 
the noninteracting Bose gas, some of the conclusions drawn here are expected to enlighten some 
aspects of the finite temperature dynamic susceptibility in weakly interacting systems, which 
can be built starting from the non-interacting case in the framework of the Random Phase 
Approximation. 
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Figure 1: Real part of xiQi^) positive u and for several values of Q. The solid, dashed, 
dotted and dash-dotted lines correspond to Q=l, 2, 3 and 5 respectively 
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Figure 2: Real part of x(q', at temperatures slightly below and above the Bose-Einstein 
transition temperature ~ 6.62. Solid line: T = 6.5, dashed line: T = 6.7, dot-dashed line: 
f = 6.8, and dot-dot-dashed line: f = 7. 
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Figure 3: Evolution of Ti?e[x(g', Y\ T)] with T and for three values of the momentum transfer q. 
The dashed, dotted, dot-dashed and dot-dashed-dashed lines stand for T = 8, 10, 50 and 500 
respectively, while the solid line corresponds to the classical limit 
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Figure 4: Momentum dependence of the qRe[x{q,Y;T)]/T at two different temperatures. The 
dashed, dot-daslied, dot-dot-dashed and dash-dash-dotted hnes stand for q = 1,5,10 and 50 
respectively, while the solid line corresponds to the limit g — ^ oo. 
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Figure 5: First moment of the real part of the dynamic susceptibihty (Eq. (p8|)) as a function of 
Q and for several values of temperature. The dotted, dashed, dot-dashed and dot-dot-dashed 
lines stand for T = 7, 8, 10 and 25, respectively, while the solid line corresponds to any T below 
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